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ON  DETERMINING  ABSOLUTE  COORDINATE  SYSTEMS  WITH  THE  AID  OF 
PLANET  OBSERVATIONS 
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BY  F.  SCHME IDLER,  MUNICH 

FROM  ASTRON.  NACHR.,  VOL.  284,  1957-59 
(RECEIVED  23  DECEMBER  1957) 

The  equations  of  condition  for  the  improvement  of  the  orbital 
elements  of  planets  from  meridian  observations  are  transformed 
in  such  a manner  that  only  the  geocentric  coordinates  ofhthe 
planets  and  of  the  sun  appear  in  the  coefficients.  For  planets, 
which  move  in  orbits  with  little  excentricity  and  inclination, 
the  equations  can  be  written  in  the  convenient  form 

An  = — E -r  b,  + a,  {t  — /,)  + a,  sin  a -i-  6,  cos  a — a,  tg  t sin  a cos  a — 6,  tg  t cos*  a 

Ad  = — D + a,  (/ — /,)  tg  c cos  a + a,  sin  a +.  (4,  + 4,  tg  *)  cos  a a,  tg  e sin  a cos  a + bi  ‘g  e cos*  « 

in  this  case,  the  coefficients  a^ , b^  are  also  a function  of 
the  quantity cc f which  is  a measure  for  the  time  period  before 
and  after  the  opposition.  The  observation  material  is  subdi- 
vided into  groups  in  wuch  a manner  that,  in  each  group,  the 
value  ofot-ix®,  and  thus  also  the  values  of  a^ , bi  are  constant. 
As  a rule,  only  the  systems  corrections  E and  D,  which  are  cha- 
racteristic for  instruments  and  observers,  are  of  interest;  the 
determination  of  the  orbital  elements  of  the  planets  can  mostly 
be  left  to  later  conversion. 

The  simplified  form  of  the  equations  makes  possible  a clear 
discussion,  as  to  which  of  the  unknowns,  under  which  conditions, 
can  be  reliably  separated.  A numerical  example  shows  that  the 
assumption  of  a small  excentricity  and  inelination  is  also 
still  permissible  for  Mars,  where  maximum  errors  should  be  ex- 
pected because  of  the  shape  of  the  orbit  and  the  nearness  of 
the  earth.  The  observation  of  large  planets  made  with  the 
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Munich  meridian  instruments  in  the  years  19^1  - 1956  are 
treated  with  the  aid  of  the  described  method. 

1.  GENERAL  COMMENTS 

The  determination  of  the  corrections  of  the  equinox  and  of  the 
equator  point,  necessary  for  the  derivation  of  absolute  star 
loci,  is  usually  carried  out  with  the  aid  of  observations  of 
the  sun  and  occasionally  also  of  the  planets  Mercury  and  Venus. 
However,  it  is  known  that  these  observations,  which  must  be 
made  during  the  day,  are  less  accurate  and  more  susceptible  to 
systematic  errors  than  observations  of  fixed  stars.  For  this 
reason,  observations  of  small  and  also  large  outer  planets  have 
been  proposed  as  alternatives  several  times  and  various  astro- 
nomical yearbooks  publish  accurate  ephemerides  of  Ceres,  Pallas, 
Juno  and  Vesta  since  1952. 

The  use  of  planet  observations  for  the  definition  of  absolute 
coordinates  is  faced  with  the  difficulty  that  no  less  than  12 
unknowns  must  be  determined.  In  addition,  the  calculation  of 
the  coefficients  of  the  equations  of  condition  is  cumbersome 
and  time-consuming.  A fictitious  example,  which  was  calculated 
by  Clmence  (l)  shows  how  extensive  the  necessary  calculation 
is.  Finally,  the  disadvantage  should  also  be  emphasized  that, 
in  the  customary  method  of  the  simultaneous  determination  of 
the  orbital  elements,  both  of  the  planet  and  of  the  earth,  as 
well  as  of  the  corrections  of  the  equinox  and  the  equator,  it 
is  not  always  possible  to  recognize  at  first  glance  which  of 
the  unknowns  are  difficult  to  separate  from  the  others. 

For  this  reason,  the  author  considered  it  desirable  to  derive 
a system  of  formulas,  which  is  easy  to  handle  numerically  and 
which  additionally  allows  a clear  overview  of  the  conditions 
of  the  separability  of  the  various  unknowns.  For  observations 
of  the  sun,  these  conditions  are  adequately  fulfilled  by  the 
known  system 


da©  =?  — E -t-  cos  e sec!d©  AL' — cos  a©  tg6©  Je  + 2 sin  a©  sec  <5©  Ah' — 2 cos  e cos  a0  sec  <5©  Ak' . 


AdQ  = — D sin  e cos  *©  J£,'  -f  sin  a©  d*  + 2 cos  a©  sin  <5©  Ah'  — 2 sin  f cos*  a©  ccs  6a  Ak 


d Ak'  j 


in  which  E and  D signify  the  sought  corrections  of  equinox  and 


equator. 


Hough  (2)  derived  a convenient  method  for  the  inner  planets 
Mercury  and  Venus,  requiring  only  the  determination  of  some 
principal  terms  of  trigonometric  series.  Numerow  (3)  discussed 
the  treatment  of  observations  of  outer  planets  in  several 
papers;  however,  the  formulas  developed  by  him  are  not  less 
cumbersome  and  vague  than  the  customary  terms  of  the  theory 
of  the  orbit  improvement.  In  contrast,  a simplification  could 
be  obtained  primarily  when  the  coefficients  of  the  equations 


of  condition  would  not  be  illustrated  as  functions  fo  the  true 


anomaly,  but  of  the  geocentric  coordinates,  because  the  latter 
can  be  taken  directly  from  the  ephemeride. 


In  the  theoretical  portion,  this  paper  has  the  objective  of 
illustrating  the  differences  between  the  observed  and  the  cal- 
culated loci  of  the  outer  planets  as  linear  functions  of  the 
unknowns  in  such  a manner  that,  if  possible,  only  the  geocentric 
right  ascension  appears  as  argument.  It  will  be  shown  that 
true  simplifications  can  be  obtained  only  by  vigorous  neglect. 

On  the  other  hand,  it  would  be  senseless  to  seek  an  accuracy 
of  more  than  10  percent;  the  differences  between  observation 
and  calculation,  with  good  ephemerides,  are  of  the  order  of 
magnitude  of  1",  and  the  average  error  of  the  unknown,  which 
are  obtained  by  the  resolution  of  the  equations  of  condition 
by  the  method  6t  the  least  squares,  are  mostly  about  io"I. 
Therefore,  it  appears  obvious  to  strike  all  terms  in  the  formu- 
las not  greater  than  10  percent  of  the  principal  terms.  In 
all  cases,  in  which  a greater  accuracy  is  required  or  is  sought 
for  specific  reasons,  it  is  probably  unavoidable  to  carry  out 
the  calculations  in  accordance  with  the  more  rigorous  formulas. 


The  neglects  which  have  been  implemented  will  not  bnly  provide 
relatively  simple  formulas,  but  it  will  also  become  clear  which 
of  the  unknowns  are  coupled  with  which  others  in  a manner  which 
is  separable  with  difficulty.  In  addition,  the  simplified 
terms  make  it  easy  to  recognize  which  of  the  unknown  can  be  de- 
termined with  centainty  only  when  the  observations  of  the 
planet  cover  a sufficiently  extensive  period  of  time  prior  to 
and  after  the  opposition.  Finally,  the  formalism  developed  in 
this  paper  provides  the  possibility  of  obtaining  at  least  as 
much  information  about  the  unknowns  from  observations,  which 
are  made  only  in  the  vicinity  of  the  opposition,  as  will  be 
possible  on  the  basis  of  the  prevailing  conditions. 

2.  INFLUENCE  OF  THE  ERRORS  OF  THE  ORBIT  ELEMENTS  OF  THE  PLANET 

Because  this  concerns  the  use  of  observations  of  equatorial 
coordinates,  all  orbit  elements  and  definitions  are  referred  to 
the  equator;  the  elements  of  the  orbit  of  the  planet  are  thus 
Mq  = average  anomaly  to  the  time  t^ 

/*■  - average  daily  motion 

2 3 2 

a = great  semiaxis,  with  pi  connected  by  /u  a = k (i+m) 

£ = sin  <p>  - excentricity 

N = ascending  node  on  the  equator 
tc  = angle  between  node  and  perihelion 
J = inclination  of  the  orbit  with  respect  to  the 
equator. 

In  the  conventional  manner,  the  heliocentric  distance  should 
finally  be  designated  as  r , the  true  anomaly  as  v and  the  geo- 
centric distance  as  0,  ; in  addition,  the  abbreviating  quantity 
u a v + w shall  be  defined  in  the  same  manner  as  in  the  eclip- 
tic calculation.  In  accordance  with  Bauschinger  (U),  the 
equations  for  the  determination  of  the  errors  of  the  orbit 
elements  are: 
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cos  dA*=  J (sin  b cos  (5  + «)  + t sin  b cos  (5  + »))  sec  <p  ( JA/#  + (<  _ /„)  4.) 

2 a i'a  r . , . 

: — jf-  y sin  *‘sin  (5  + u)  4u 

+ — (sin  £ sin  6 cos  (£  -f  u)  — cos  <p  sin  b sin  (B  -f  v))  dip 
+ j sin  b cos  (B  + u)  ds  + ~ cos  b sin  u dj  — ^ cos  b cos  u sin  J dN  , 

a 

Ah  = - (sin  c cos  (c  + ")  + <•  sin  f cos  (C  + if1))  sec  <p  (dM0  + (l  — 1„)  d/s) 

2 a la  r . 

3*  '?Sln  C Sm  w)  Ap 

+ ^ (sin  E sin  <•"  cos  (r  + u)  - COS  <p  sin  c sin  (C  + u ))  d<f 

f f y • • 

+ J sin  c cos  (C  + u)  da  + — cos  c sin  u dj  — — cos  c cos  u sin  J dN . 


In  this  case,  instead  of  the  quantities  Jp,  Aq,  used  by  Bau- 
schinger,  the  direct  corrections  dN  of  the  elements  are 

introduced;  the  following  relationships  apply 
dp  - : sin  J sin  w dN  -i  co<  \]  . 

dg  -=i  sin  ] cos  »•'  AX  — >in  i.  \]  , ^ 

J*  = cos  / dN  ~ tlw . 


The  auxiliary  quantities  b,  B are  defined  by  the  system 


sin  b sin  B — — sin  (n  — X)  , 
sin  b cos  li  = cos  / cos  (x  — .V) , 
cos  6 — — sin  J cos  (s  — N) 


(4) 


while,  for  c,  C,  the  relationships 

Mil  r sin  C - — sill  A co>  (s  — .V)  . 
sin  c co-  C ■ sin  J cos  A — cos  J sin  A sin  ( \ - .\ ) , 

cos  c - - cos  ] cos  A -t  sin  J sin  A sin  (*  — N) 


apply  as  definition. 


In  order  to  avoid  the  known  uncertainty  for  small  excentrici- 
ties  and  inclinations,  the  new  unknowns  are  better  used 


t 


y 

i 


dLt  = sec’  <p  d\f0  + ds 
dv  — c ds  , 
dK  = sin  J dN  . 


(6) 


In  (2),  the  semimajor  axis  a ini  the  coefficients  of 
should  additionally  be  replaced  by  the  expression  r 1 yt  (i ♦ ^ 

and,-  in  the  coefficients  of  A <p  the  factor  a sin  E by  the  ex- 
pression v if  v . If  u - v is  finally  also  written  in  (2) 
instead  of  w,  the  following  form  of  the  equations  of  condition 
is  obtained: 


cos  v)  («  sin  v /JL  - Sin  v Ay)--  cos  r cos  v A<p  . 


p = H J/4  + (l  + * 
3* 

f)  = AL  + (2  + e cos  v)  (< ' 


, AT  — cos  v /Jvi  4-  (sec 


/-/-vc  rr  I «iin  V 


(7) 


AL  = AL,+  (<  — gsec*?>^- 


In  equations  (7),  r,  u and  v must  be  expressed  by  the  geocen- 
tric coordinates  of  the  planet  and  of  the  sun.  It  is  known 
that 


x = o cos  a cos  <5  — 7?  cos  «0  cos  <30  = r (cos  AT  cos  u — sin  AT  sin  u cos  /) , 
y = 0 sin  <x  cos  <5  — 7?  sin  ctQ  cos  J0  = >-  (sin  N cos  u + cos  N sin  u cos  J) , 
z = 0 sin  6 — R sin  60  = r sin  u sin  J , 


where  R,  J & are  the  geocentric  coordinates  of  the  sun. 

If  the  two  auxiliary  quantities 


R cos  <5q 
0 cos  6 


cos  (a  — a0) , 


V = 


R cos  <50 
q cos  d 


sin  (a  — a. 3) 


(9) 


are  introduced,  these  relations  are  obtained  from  (8) 


r sin  u = 0 (1  — |)  cos  5 sin  (a  — A7)  sec  J 4-  0 jj  cos  <5  cos  (a  — N)  sec  J , 
r cos  w = g (1  — f)  cos  <5  cos  (a  — -V)  — 0 rj  cos  <5  sin  (a  — AT)  • 


(10) 


Formulas  (10)  make  possible  the  calculation  of  r,  u and  v = 
u - w from  the  geocentric  coordinates  of  the  planet,  which  can 
be  taken  from  the  ephemeride.  The  two  parameters  £ and  are 
essentially  a function  of  ot-x.@  and  are  a measure  of  the  time 
interval  of  the  observations  from  the  opposition. 


Because  most  planets  move  reasonably  close  to  the  ecliptic,  it 
is  useful  to  express  sin  S as  a function  of  x . Using  the 
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second  equation  (10),  the  following  relationship  is  found  from 
the  third  equation  (8) 

e sin  b = o (i  - £)  tg  / cos  b sin  (a  - -V)  + 2 >>  tg  / cos  A cos  (*  - M -r  * *™  * s • 

Because  the  known  relationship  exists,  the 

following  is  obtained 

/?  sin  = R tg  t cos  (5©  sin  a©  = g tg  f cos  b (c  sin  a — r\  cos  a) . 

If,  herein,  «tis  replaced  by  N + (m.  - N),  the  ft  .lowing  expres- 
sion is  obtained  for  sin  <T 

sin  b --  tg  ] cos  b sin  (a  — .V)  — " 

with  I (1I) 

;•  =-*  (tg  J — cos  .V  tg  t)  — »/  sin  .V  tg  r)  cos  6 sin  (a  — A ) — (//  (tg  ] — cos  A tg  fi 
-r  ; sin  .V  tg  t.)  cos  b cos  (a  — .V) . _ 

For  planets,  which  move  exactly  in  the  ecliptic,  J * i and  N = 0, 
and  consequently  also  = 0.  For  planets,  of  the  which  the 
orbit  inclination  with  respect  to  the  ecliptic  is  small,  y can 
be  considered  as  a small  quantity  and  can  be  neglected  complete- 
ly when  the  requirements  for  accuracy  are  small. 

The  advantage  of  relationship  (ll)  for  sin  <f  consists  of  the 
fact  it  allows  the  expression  of  the  quantities  c,  C as  func- 
tions of  b,  B and  do  therefore  not  have  to  be  calculated  sepa- 
rately. If  expression  (ll)  for  sin  J is  inserted  in  the  two 
first  equations  (5),  the  following  is  found 

.-in  <•  sin  ('.  - ig  ] cos  b cos  (a  — A")  sin  h >in  If  ;•  sec  J sin  h cos  Ji , 

'in  c cos  C = tg  ] cos  <)  cos  (a  — A’)  sin  b Co*  If  — ;•  cos  ] sin  h sin  11 . 

In  the  terms,  which  are  multiplied  with  the  small  quantity 
sec  J - cos  J = 1 can  be  inserted,  and  it  follows 

r sin  c sin  (C  - »<)  = r tg  ] cos  b cos  (a  — .V)  sin  b sin  (13  -f  u)  + yr  sin  b cos  (B  -f  «) . 1 £: 

r sin  c cos  (C  - «)  = r tg  / cos  b cos  (a  — .V)  sin  b cos  (B  -f  u)  — r sin  b sin  (B  + u) . J 

In  addition,  an  expression  for  cos  c can  be  found,  when  ( 11 ) is 
inserted  in  the  third  equation  (5) 

cos  c = cos  / cos  b -j-  sin  / tg  / cos  b sin5  (a  — A’)  — •/  sin  J sin  (a  A') . 

By  neglecting  the  square  of  ^ , the  following  can  be  written 
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l 


cos  c = cos  / cos  6(1  + tg*  6)  + y sin  / sin  (<x-N)=  cos  / sec  <5  + ;<  sin / sin  (a  — JV) . ^ 

All  quantities,  which  occur  in  the  equations  of  condition  (7)> 

can  be  calculated  with  the  aid  of  formulas  ( h ) , (9)*  (10),  (12) 

and  (13).  The  derived  formulas  are  completely  formal,  except 

that,  in  the  terms  multiplied  with  the  small  quantity.)',  sec  J 

2 

= cos  J = 1 is  inserted  and  ^ is  neglected.  In  most  cases, 
jj  = 0 can  also  be  assumed. 


3.  INFLUENCE  OF  THE  ERRORS  OF  THE  ORBIT  ELEMENTS  OF  THE  EARTH 

As  a rule,  only  four  earth  orbit  elements  are  considered  as 
requiring  improvement.  The  following  unknowns  are  used 


AL* 

= correction 

of 

the 

A c 

= correction 

of 

the 

Ah' 

= correction 

of 

the 

A*' 

» correction 

of 

the 

mean  longitude  of  the  sun 
skew  of  the  ecliptic 
quantity  h'  = e'  cos  nr' 
quantity  k'  = e'  sin  n 1 , 


Newcomb  (5)  derived  the  formulas,  which  represent  the  influence 
of  the  corrections  of  the  earth  orbit  elements  on  the  planet 
loci.  However,  an  independent  treatment  of  the  problem  appears 
to  be  more  suitable  for  the  purposes  of  the  present  work.  As 
a starting  point,  equations  (8)  can  be  used  in  the  form 


p cos  a cos  6 = x + x' , 
g sin  a cos  6 = y 4-  y' , 
g sin  6 = z + i' 


in  which  x',  y',  z'  are  the  geocentril  rectangular  equatorial 
coordinates  of  the  sun.  They  are  given  by  the  expressions 


x'  — R cos  aQ  cos  dQ  = R cos  A©  • 

y'  = R sin  a©  cos  6q  = R cos  e sin  A©  , 

z’  = R sin  3©  — R sin  e sin  A©  . 


(15) 


of  Aq t the  following  differential  equations  are  obtained  in  a 
known  manner  from  equations  (lU) 

o cos  6 da  =•  cos  a dy'  — sin  ft  Jx' , 

q AS  = cos  <5  A:’  — cos  a sin  6 Ax'  — sin  a sin  S Ay'. 

However,  since  the  relation  z'  = y'  tanf  applies  in  accordance 
with  equations  (15)  > A2'  can  be  replaced  by Aj 1 
and  the  following  is  found 


n cos  J .1a--=  cos  a dy*  — sin  a Jx' , 

n .16  =■  v’  socs  f cos  i)  A*  — cos  -v  sin  6 l v'  — (ty  r cos 


’I 

:os  6 — s;n  i sin  6)  Jy' . | 


With  the  aid  of  equations  (15),  Ax'  and  Ay'  can  be  expressed 
as  functions  of  A R , A<^  o and  > and 

.-IA’ 

1'  — ' ,,  — y see  c l/.  j , 

4 . ..IK  , , ,,  , (*7) 

dy  - y w 4-  X COS  f JA-  — y'  tg  c . 
is  obtained. 

The  interrelation  between  the  coordinates  R,  A <£>  and  the  orbit 
elements  of  the  earth  is  given  by  the  known  formulas 

A’  «'  (i  --  /i’  cos/.j  - - //  >in/  .)  . /„  - ; /»'  sinAw  — 2 /.•'  cos 

in  which  the  square  of  the  excentricity  of  the  earth  orbit  is 
neglected.  Through  differentiation 

A li  — x’  Ah’  — y’  sec  t \k' , | 

J/i,  --  At'  -(-  see t Ah'  — Ah'.  I ^ 

is  obtained. 

These  expressions  must  be  inserted  in  (17)  and  the  result  is 


Ax’  « — y'  see  e AL'-X-  j/,-  + J(t>( 

**  A 

dy'  - x'  cos  £ AL’  + Ah’  - tj^.L  + If'!  c°li  Ak'  _v> 


Ak'  — y'  tg  e Ac. 


These  expressions  must  be  inserted  in  (l6)  in  order  to  obtain 
the  illustration  of  Ax  and  as  functions  of  the  improvements 
of  the  earth  orbit  elements.  The  result  of  the  substitution 
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shall  not  be  calculated  in  greater  detail  at  this  point,  because 
the  formulas  become  very  cumbersome,  but  they  do  take  on  a 
simpler  form  with  the  aid  of  the  neglect  in  the  later  part  of 
this  paper.  However,  it  should  be  noted  that  the  coefficients 
in  (19)  are  fundamentally  a function  of  the  rectangular  coordi- 
nates of  the  sun,  which  can  be  taken  from  yearbooks. 

U.  INTRODUCTION  OF  SIMPLIFICATIONS 

The  formulas  developed  under  2.  and  3-  provide  the  coefficients 

of  the  equations  of  condition  as  functions  of  the  geocentric 

coordinates  of  the  planet  and  of  the  sun;  the  calculation  of 

the  locus  of  the  planet  in  its  orbit  can  be  avoided  with  their 

aid.  Yet,  the  required  calculations  should  not  be  significantly 

less  than  with  the  use  of  the  conventional  system  (2).  However, 

a considerable  simplification  can  be  obtained  if  the  excentri- 

cities  and  the  square  of  the  inclination  are  neglected.  Because 

most  planets  move  in  the  vicinity  of  the  ecliptic,  J is  approxi- 

2 

mately  £ and  the  assumption  J = 0 is  just  tolerable  with  the 
requirement  for  an  accuracy  of  about  10  percent,  which  was  dis- 
cussed in  the  introduction. 

2 

When  J is  neglected,  sec  J = cos  J = 1 and  the  following  is 
found  from  equations  (U)  and  (10) 

r sin  b sin  (Z?  4-  i<)  = or,  cos  6 , 
r sin  b cos  (B  -J-  «)  =5(1  — f)  cos  d . 

Because  of  J ~£  , the  difference  J - £ is  of  the  order  of  magni- 
2 

tude  J and  can  be  neglected;  for  the  same  reason,  N as  0,  and 
the  product  sin  N tan  £ can  be  neglected  as  a quantity  of  the 
second  order.  Finally,  under  the  assumed  prerequisites,  the 
declinations  of  the  sun  and  of  the  planet  are  greater  than  £ 
by  at  most  insignificant  amounts  and,  therefore,  cos  S = cos  a 
= 1 can  be  used.  Under  consideration  of  (20),  the  following 
relations  are  then  obtained  from  ( 1 1 ) and  (12) 
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sin  6 = tg  e sin  a . 

r sin  c sin  (C  -f  «)  = r tg  e cos  a sin  b sin  (B  -f  11)  = o t]  tg  s cos  a , ) (21) 

r sin  c cos  (C  4-  11)  = r tg  e cos  a sin  b cos  (B  4-  n)  = o (1  — f)  tg  £ cos  a . j 

With  the  aid  of  the  same  neglects,  the  expressions 

r cos  b sin  11 
r cos  b cos  11 
r cos  c sin  » 
r cos  c cos  11 

can  be  derived  from  (4),  (10)  and  (13),  in  which  the  auxiliary 
quant  it  ies^'have  the  significance 

Pi  — (1  — c)  cos  .V  4-  ij  sin  A’ , | 

— — (t  — i)  sin  .V  4-  > ] cos  .V  J (_^) 

With  the  assumptions  r = a,  sec  J = cos  S = 1 and  under  consi- 
deration of  v = u - w,  the  following  expressions 

sin  ;•  = rs  sin  a — r cos  a , 
cos  v — t sin  \ 4-  a cos  a , 

are  obtained  from  equations  (10)  for  the  true  anomaly  v,  which 
occurs  in  the  equations  of  condition  (7),  in  which  the  auxiliary 
quantities  6 and  z are  defined  by 


= — 0 tg  e cos  a (pi  sin  a 4-  ?t  cos  a)  , 
= — n tg  £ cos  a (pi  cos  a — sin  a)  , 
= o (pi  sin  a 4-  9i  cos  a)  , 

= n (/>,  cos  a — <h  sin  a) 


n = -J  (t1  — .*)  cos  (-•  -r  .V)  4-  r,  Sin  (;r  4-  .V)) , 
r = ~ ((r  — f)  sin  (v  4-  A’)  — >j  cos  (w_  4-  A*)) . 


(-’4) 


If  expressions  (20),  (21),  (22)  and  (24),  which  are  obtained  by 
simplifications,  are  inserted  in  the  equations  of  condition  (7), 
then  the  following  illustration  for  Aon  and results 

Aa  — At(l  — l0)  4-  B0  4-  -4,  sin  a 4-  Bj  cos  a — A . tg  e sin  a cos  a — B.  tg  e cos*  a , 1 
Ad  = (,40  (/ — /„)  4-  B#  4-  Al  sin  a 4*  B1  cos  a)  tg  £ cos  a 4-  At  sin  a 4*  B.  cos  * • J 


In  this  case,  the  coefficients  A ^ , B are  related  to  the  sought 
corrections  of  the  orbit  elements  by  the  following  expressions 


I 


a 


At=  ( i — f)  sec*  <P  Afi . \ Bt  = (i  — f)  JL0  — — r)  A/i , 

At  = —r)  (o  Av  + t A<r)  — 2 {1  —i)  (t  Av—oA<p),  j B,  = 7](rAv— a Ay)  — 2(1—  {)  (a  Av  + rAp). 
A*  = PiAJ  + ?!  AK . Bt  = qxAJ  — piAK. 


(21 


All  auxiliary  quantities  are  defined  by  equations  (9).  (23) 
and  ( 2k ) . 


Equations  (16)  and  (19),  which  describe  the  influence  of  the 
errors  of  the  orbit  elements  of  the  earth,  can  be  simplified 
by  the  same  neglects.  Under  the  assumption  of  sec/  = cos  £ 
= 1 and  using  the  expression  ( 2 1 ) for  sin 
o Ax  = cos  a ,dy'  — sin  a Ax' , 

q Ad  = y'  At  — tg  £ cos  a (sin  a Ax'  — cos  a Ay')  — y’  Ae  + Q tg  4.  cos  a Ax  . 
is  obtained  from  (l6).  Expressions 

x'  = 0 cos  <5  (£  cos  a + r,  sin  a) . y'  — e cos  <5  (f  s‘n  ® *7  cos  *)  • (2®) 

are  found  from  (9)  and  (15)  for  the  rectangular  coordinates 
of  the  sun  x ' , y ' . 


r 2 

If  now  again  cos  = cos  c = 1 and  tan  c = 0 , the  following 
illustration  of  the  deviations  caused  by  the  errors  of  the 
earth  orbit  elements,  which  should  now  be  designated  as  &*.' 
and  AS  ' , results  from  equations  (19),  (27)  and  (28): 


Ax'  = £AL'  + +R2(y'  Ah'  + ’7--~  -R2  ( X'  Ah’  — (f  sin  a — rj  cos  a)  tg  t cos  a As',  j 

Ad’  = (i  AL'  -f  - — ^ ^ Ah'  4-  — — r ^ ^ *6  6 cos  <*  + (f  sin  « — t]  cos  a)  Ae . 


(29) 


Using  the  relations  (28),  the  coefficients  of  these  equations 
can  be  expressed  as  functions  of  ^ , as  well  as  of  £ and  n 


rjx'  + 2Sy' 
R 

T)y'  — 2(x' 
R 


(2  f*  + rj1)  sin  a — ^ (r)  cos  a ; 


R 

2 

R 


— -|  (2  + *7*)  cos  * — £ 7 sin  a . 


It  can  be  seen  from  equations  (29)  and  (30)  that  the  deviations 
Aa'  and  , which  are  a consequence  of  the  errors  of  the  orbit 
elements  of  the  earth,  can  be  illustrated  in  the  same  general 
form  (25)  as  the  quantities  4*and  Af , which  are  based  on  the 
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errors  of  the  orbit  elements  of  the  planet.  If,  additionally, 
the  corrections  E and  D for  equinox  and  equator  are  still  added, 
the  result  is  the  general  illustration  for  the  difference 
between  the  observed  and  the  calculated  right  ascensions  and 
declinations  of  a planet 

Aa  = — E + 6,  + a,  [t  — tt)  + a,  sin  a ■+■  6,  cos  a — a,  tg  4 sin  a cos  a — 6.  tg  t cos  a* , I ^i) 

Ad  — — D + a#(/ — /,)  tg  ecosa  -fa,  sin  a + (6,  -f  60tgf)  cosa  -f  tg  f sin  a cos  a + ijtgtcos'a.  J 


The  significance  of  the  coefficients  a^  , b^  follows  from  equa- 
tions (26),  (29)  and  (30)  through  the  addition  of  the  quanti- 
ties derived  from  the  errors  of  the  orbits  of  the  planet,  res- 
pectively of  the  earth 


at  = [J—  f)  SPC*  7 Aft  , 

K = (1  - £)  ALt  — ™15 ry  A,,  4-  £ AV , 

а,  = — jy  (a  Av  J-  t Jf/j  — 2 (1  — £)  (t  Ar  — <r  .1?)  J-  £ (2  « -f  »/=)  Ah'  — £ r)  Ah'  . 

б,  = »y  (t  Av  — o A<f)  — 2 (1  — £)  (<r  /J>-  + t Jy)  — £ £ »y  JV  — ^ (>  « -f  iy3)  JAr' . 

a,  = J/  + y,  AK  -f  £ Js  , 

b.  = y,  d/  — />,  J A"  — t]  Ac . 


If  it  is  additionally  desired  to  account  for  a correction  A p 
of  the  precision  constant  and  a time  dependance  of  AL’  in  the 
form  AL’  = Aq  + /'^(t  - t^),  then  the  expressions  Av  + 
must  be  added  at  the  right  in  the  first  equation  (32)  and  ^L1 
must  be  replaced  by  A q in  the  second. 

The  advantage  of  equations  (31)  consists  of  the  fact  that  the 
coefficients  a^  , b,.  are  constant  when  ct-  - has  the  same 
value  for  all  observations.  Actually,  with  the  neglects  which 
are  introduced,  the  quantities  <?  , f andij  are  only  a function  of 
the  semimajor  axes  of  the  two  orbits  and  of  - *-<3 . With  the 
assumptions  cos  J = cos  (5"  = cos  <T  43  * 1 and  e = e 1 =0,  conse- 
quently r = a,  the  relations 

y3  — 2 R q cos  (.t  — «©)  -j-  Rz  = a1 , 

t R . R . (33) 

£ = — cos  (a  — ««) , n = — sin  (.t  — a©) . 

e " e 

are  obtained  from  equations  (9)  and  (10). 


bzkazjX, i 1 .. * - 
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Therefore,  for  all  observations,  for  which  the  value  X - <*  a is 
equal  or  falls  into  a sufficiently  small  interval,  the  coeffi- 
cients a.^  , b^  must  be  determined  through  the  solution  of  the 
system  (31)*  The  unknowns  themselves  then  follow  from  (32). 

In  this  manner,  the  determination  of  the  12  unknowns  is  sub- 
divided into  several  small  steps. 

5.  THE  DISCUSSION  OF  THE  FORMULAS 

Because  the  meridian  observer  is  mostly  only  interested  in  the 
systems  improvements  E and  D,  the  separability  of  these  two 
quantities  shall  be  discussed  first.  The  determination  of  D 
is  simplest;  it  is  equal  to  the  constant  term  in  the  equation 
for  £if ; uncertainties  cannot  arise  when  the  material  of  the 
observations  is  good  and  quite  well  distributed.  A division 
of  the  material  in  accordance  with  the  value  is  not  even 

absolutely  necessary;  fitting  all  observations  in  accordance 
with  the  second  equation  ( 31 ) , values  are  obtained  for  the 
quantities  a^ , b^ , which  represent  an  approximate  average. 

Again  assuming  a reasonably  uniform  distribution  of  the  entire 
material  of  the  observations,  the  numerical  value  of  D will 
be  adulterated  by  insignificant  amounts  at  most. 

Contrary  to  the  occasionally  expressed  opinion,  it  must  be 
emphasized  that  the  equator  correction  is  not  simply  equal  to 
the  arithmetic  mean  of  all  Because  of  the  factor  cos‘  <*  , 

the  last  terms  in  (31)  always  have  the  same  sign  and  definite- 
ly contribute  to  the  arithmetic  mean  of  all  observed  values. 

The  determination  of  the  equinox  correction  E,  which  occurs  in 
(31)  as  the  fixed  relationship  -E+b^,  presents  somewhat  greater 
difficulties  in  principle.  The  quantity  bQ  must  therefore  be 
determined  from  the  declinations,  in  which  it  occurs  multi- 
plied by  tan  £,  however.  For  this  reason,  E,  in  comparison 
with  the  other  unknowns,  is  encumbered  with  an  uncertainty, 
which  is  increased  by  the  factor  cotan  £'.  These  subjects  are 
not  new,  but  they  can  be  seen  particularly  clearly  with  the 


mmmm 


simple  formulas  (31) 
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In  addition,  the  determination  of  E poses  high  requirements 

on  the  quality  of  the  observed  right  ascensions,  from  which 

the  reliable  value  of  b2  must  be  determined;  the  term  with 

bQtanfcan  only  then  be  separated  from  bg  from  the  declinations. 

Because  b?  is  associated  with  the  mostly  small  factor  tan  £ 

2 d 

cos  i*  in  the  right  ascensions,  an  additional  uncertainty 
prevails  in  this  manner. 

Although  the  material  on  the  observed  right  ascensions  raigght 
not  fulfill  the  conditions  for  a reliable  determination  of  b2, 
but  has  a sufficiently  good  distribution  in  the  variable 

the  dependence  of  the  coefficients  ai  , b^  from  this 
second  variable  can  be  utilized  in  order  to  obtain  a reliable 
value  of  E.  In  this  case,  the  three  quantities  a2,  b^  = 
bg  + bg  tanf  and  EQ  = -E+bQ  can  be  determined  with  the  full 
reliability,  which  is  associated  with  the  observations.  From 
( 32),  the  relationships 

Pi  AJ  -j-  ?!  AK  -f  $ Ae  = a,  , 

?i  AJ  — />!  JA'  — t]  An  -I-  E tg  £ = 6,  — E„  tg  e 

can  be  derived  and  the  four  unknowns  4j,  AK , and  E can  be 
determined  when  the  material  allows  a subdivision  into  at  least 
two  considerably  different  values  of  *--*.0.  This  possibility 
was  largely  used  in  the  evaluation  of  the  Munich  planetary 
observations,  which  are  further  reported  below. 

The  determination  of  the  equinox  correction  E is  configured 
still  more  advantageous,  when  the  errors  of  the  orbit  elements 
of  the  earth  are  either  neglected  or  can  otherwise  be  assumed 
as  known.  In  accordance  with  the  second  equation  (32),  the 
simple  relationship  for  the  quantity  -E+bQ,  which  is  determin- 
able from  the  right  ascensions  alone 

E ~ b0  = Ea  = E — (1  i)  ALg  ^ - i/  Aft , (35) 

3 « 

from  which  E can  be  determined,  if  the  material  of  the 
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observations  is  distributed  over  at  least  two  significantly 
different  values  of  and  Vis  taken  from  the  first 

equation  ( 32 ) . 

The  corrections  of  the  orbit  elements  of  the  earth  and  of  the 
planet  can  be  found  from  (32).  Again,  as  can  be  seen  immedia- 
tely, the  numerical  values  a^ , b^  are  required  for  at  least 
two  significantly  different  values  - * o!  -he  occupation  of 
as  wide  a range  in  this  variable  is  of  course  desirable.  If 
the  corrections  of  the  earth  orbit  elements  can  be  neglected, 
the  problem  becomes  correspondingly  simpler. 

The  simplification,  which  can  be  obtained  through  the  neglect 
of  the  corrections  of  the  earth  orbit,  has  considerable  signi- 
ficance, primarily  for  the  outer  planets,  which  are  at  a con- 
siderable distance.  It  is  obvious  that  the  observations  of 
these  planets  can  provide  little  information  about  the  orbit 
elements  of  the  earth;  analytically,  this  is  expressed  by  the 
fact  that,  in  (32),  all  corrections  of  the  earth  orbit  are 
multiplied  with  the  factors  £ or  A , which,  corresponding  to 
expressions  (33),  become  very  small  for  very  distant  planets. 

Finally,  it  is  possible  to  continued  in  so  far  as  the  systema- 
tic errors  of  the  instrument  can  be  considered  as  periodic 
functions  of  the  right  ascensions,  thus  the  errors  A and 
/3c*  can  be  introduced  in  the  form  of  terms  which  are  propor- 
tional to  sinX  and  cos*  . The  amplitude  of  these  terms  would 
have  be  added  in  the  system  (32)  in  the  appropriate  equations. 

It  can  be  seen  that  these  terms  can  be  reliably  separated  from 
the  other  unknown  only  when  the  observations  uniformly  cover 
a sufficiently  broad  range  of  the  variables  <*.  - *.  0 ; this  is 
equal  to  the  requirement  that  the  observations  must  extend  over 
a sufficiently  long  time  period  prior  to  and  after  the  opposition. 
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6.  VERIFICATION  OF  THE  ADMISSIBILITY  OF  THE  SIMPLIFICATIONS 


It  was  tested  on  the  basis  of  a numerical  example  whether  the 
neglect  of  the  excentricities  and  the  squares  of  the  inclina- 
tions is  permissible.  As  the  specimen  calculation,  the  obser- 
vations of  Mars,  made  in  the  years  19l*l  to  1956  with  the  Munich 
meridian  instruments,  were  selected,  because  the  greatest  errors 
are  suspected  in  the  case  of  this  planet  because  of  its  excen- 
tricity  and  nearness  to  the  earth.  All  declinations  are  pub- 
lished in  the  catalog  of  the  author  (6)  on  azimuth  circle  ob- 
servations; however,  the  values  which  are  given  there  are 
already  corrected  because  of  the  systems  corrections,  which  are 
here  investigated  and  are  therefore  different  from  the  values 
given  in  this  paper.  Up  to  the  beginning  of  1950,  the  right 
ascensions  are  published  by  Labitzke  (7),  (8);  in  addition, 
the  author  owes  thanks  to  Mr.  Labitzke  for  the  premature  re- 
lease of  unpublished  observations  of  the  years  1950-1952. 

Table  1 shows  the  normal  loci,  which  are  formed  from  the  ma- 
terial. Three  groups  were  formed,  one  each  for  observations 
with  small  (A),  medium  (B)  and  large  (C)  value  of  - <*  c s. 

In  the  column  "location",  which  occurs  only  in  the  declinations, 
M = Munich  and  C = Canberra  because,  in  195^/55,  observations 
were  made  in  Canberra,  Australia,  using  the  Munich  azimuth 
circle.  All  other  designations  in  Table  1 do  not  require  an 
explanation. 

In  the  distribution  of  the  observations  into  groups,  no  fixed 
ranges  of  - X g were  maintained,  but  individual  observations 
were  absorbed  in  another  group  for  the*  purpose  of  a more  advan- 
tageous distribution.  It  proved  to  be  impossible  to  form  more 
than  three  groups,  which  entailed  considerable  advantages. 

The  improvements  of  the  orbit  elements  of  the  earth  can  be  ig- 
nored for  the  questions  which  here  alone  are  discussed , whether 
the  leglects  which  have  been  carried  out  are  permissible.  The 
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Table  1.  Normal  Loci  of  Mars 
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equinox  corrections  E were  also  neglected,  so  that,  in  addition 
to  the  six  improvements  of  the  orbit  elements,  two  additional 
unknowns  remained,  namely  the  equator  corrections  for  the  azi- 
muth circle  observations  in  Munich  and  Canberra,  which  can  of 
course  be  different  and  were  different.  Eight  unknowns  must 
therefore  be  determined,  which  were  determined,  on  the  one  hand 
by  the  method  of  approximation  described  in  this  paper  and,  in 
addition,  by  the  application  of  the  formal  system  (2)  from  the 
normal  loci  of  Table  1.  As  the  mean  value  of  * - ^ , £ 

and  ^ the  following  resulted  for  the  three  groups 

.Group  A a— = 7h36"\  n ~ o.8S,  £ = — o.qf),  j;  = +1.04 

5 0.69  —1.05  +T.00 

5 0.55  —1.77  -fo.43 


The  determination  of  a^ , Ik  was  carried  out  by  way  of  approxi- 
mations. The  second  equation  (31)  was  written  in  the  form 

a,  sin  a + b3  cos  a — D — dd  — aa  ( t — 19)  tg  e cos  « — at  tg  e sin  a.  cos  a — bx  tg  e cos*  a.  (36) 

and  was  used  for  the  determination  of  the  four  unknowns  a.^ , 
b^,  Dm  and  Dc , under  the  hypothesis  for  the  quantities  a^  , a^ 
and  b1  (in  the  first  approximation,  these  quantities  were  set 
equal  to  zero),  for  the  observations  of  each  group.  With  the 
values  which  were  thus  obtained,  the  system 

E*  + *0(t  — <#)  4-  «isina  + ^cosa  — 6jtgecos*«  = Ao.  + a,  tg  e sm  a cos  a , 1 (37) 

«o  ( t — tt)  tg  e cos  a -f  ax  tg  c sin  a cos  a + tg  e cos*  a = A&  + D a,  sin  a j cos  a J 

was  solved  for  the  five  unknowns  EQ,  a^,  a^ , b^  and  bg , in  that, 
on  the  right  sides,  the  correction  terms  were  calculated  with 
the  results  of  the  previous  approximation.  If  the  material  of 
the  observations  on  the  right  ascensions  were  quantitatively 
equal  to  the  declinations,  it  would  be  better  to  use  only  the 
first  equation  (37)  for  the  determination  of  the  five  unknowns; 
however,  it  proved  to  be  necessary  to  also  utilize  the  second 
equation.  The  method  was  iterated  until  the  values  of  the 
unknowns  no  longer  changed. 

In  the  determination  of  the  two  quantities  D from  (36),  it  must 
be  observed  that,  while  the  coefficients  and  b^  are  differ- 
ent in  each  group,  the  values  of  D,.  and  l>-  must  however  be 
equal.  In  order  to  obtain  the  best  value  of  the  two  D within 
the  scope  of  the  error  theory,  the  normal  equations  which  were 
obtained  from  (36),  were  reduced  in  each  group  by  elimination 
until  only  the  two  unknowns  and  remained;  these  equations 
from  all  three  groups  were  added  and  their  solution  resulted 
in  the  associated  values  of  and  Dc  in  each  stage  of  the 
approximation . 

The  formal  calculation  in  accordance  with  formulas  (2)  re- 
sulted in  the  following  system  of  equations  of  condition: 


41.41  d.v0 

— 0.81  d// 
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— 0.3b  d/ 
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Dm 

— 0.65 

— 0.07 

—0.76 
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— 
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— 0.61 
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— 

Dm 

— 0.5S 

0.00 
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-7- 1.7" 
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—0.65 
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-1  53 
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— 
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4 1.69 

4 2.0S 

40.51 

-139 

—0.51 

—324  = 

—0.55 
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Dm 

4-O.Q2 

40.77 

40.77 

4-0.7S 

40.69 

— 2.02  = 

— 0.20 

— 

Dm 

40.54 

-0.3S 

4o.S(i 

--0.4S 

41.00 
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40.50 

(3») 


The  result  showed,  from  the  approximate  calculation,  that  the 
two  equator  corrections  for  the  azimuth  circle  observations 
had  the  values  DM  = +0"02  and  Dc  = +0"31.  Of  the  other  unknowns, 
^2  very  unreliable  in  all  cases  and  could  not  be  determined 

at  all  for  groups  3 and  C.  The  results  are 

aa  ai  at  6,  4,  l3 

GroppeA  +o‘.'ii  4o"24  — 2,-42  +o"6i  +177  -f-o"68:  — o'bs 

B 3.06  4o.o8  +0.S3  — 0.01  -i-0.53  unb.  — 0.45 

C —1.95  —0.06  — 0.92  — 0.56  —0.34  unb.  —0.55 

Neglecting  the  corrections  of  the  earth  orbit  and  of  E,  the 

improvements  in  the  orbit  elements  of  Mars  must  be  found,  in 
accordance  with  (32)  from  these  values  of  a^ , b..  In  Table  2, 
the  results  are  compared  with  those  of  the  solution  of  the 
completely  formal  system  (38). 
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Table  2.  Comparison  of  Approximated  and  Formal  Results 


(j)  Unbckannte 

strenger^ 

Wert 

genihertet^ 
Wert  (X 

AL, 

a a 

-o  57±°  94 

a a 

—o. 82  + 0. 66 

A/i 

— 007±0.24 

+O.OI  ±0.23 

Jr 

—0.30±0.12 

— I o.  >5*0.19 

A<p 

— 0.2710.14 

— o 0510.19 

+0.11  + 0.16 

+o.o3±o.i4 

AK 

o.oo±o.i8 

— < 0.1710.13 

D„ 

— O.XO±O.I9 

+0.0210.14 

DC 

+0.I2±0.52 

+0.3110.40 

Key : 

1.  Unknown 

2.  Formal  value  * 

3.  Approximated  value 

With  consideration  of  the  mean  error,  the  agreement  is  good 
in  all  cases.  Even  if,  with  three  times  the  observation  ma- 
terial, the  mean  error  would  be  smaller  by  a corresponding  fac- 
tor, factual  differences  would  still  not  occur  between  the  two 
solutions . 

The  fact  that,  in  the  case  of  most  unknowns,  the  approximation 
method  results  in  considerably  smaller  mean  errors  is  of  course 
only  apparent.  The  approximation  method  is  based  on  the  solu- 
tion of  systems  with  few  unknowns,  where,  formally,  the  mean 
errors  must  be  small.  It  is  however  known  that,  as  a rule, 
these  methods  provide  fairly  correct  values  of  the  unknowns, 
even  though  their  accuracy  is  stated  too  positively.  In  con- 
clusion, it  can  be  said  that  the  neglects,  on  which  the  approxi- 
mation methods  is  based,  appear  Justified. 

7.  EVALUATION  OF  THE  MUNICH  OBSERVATIONS  OF  LARGE  PLANETS 

The  observations  of  Jupiter  and  Saturn,  which  were  made  with 
the  Munich  instruments,  were  also  evaluated  in  accordance  with 
the  described  method.  Because  of  the  short  orbital  arc,  which 
was  passed  through  in  the  years  19**1  to  1956,  only  a summary 
treatment  was  possible  for  Uranus  and  Neptune.  In  the  subse- 
quent work,  Mr.  Petri  reports  on  the  use  of  the  observations 
of  the  small  planets  Ceres,  Pallas,  Juno  and  Vesta. 


The  observations  are  published  in  the  same  publications  (6), 
(7),-  (8)  as  the  Mars  observations  referred  to  above.  The  de- 
clinations observed  on  the  meridian  circle  were  not  used  be- 
cause of  the  accuracy  which  is  inferior  to  that  of  the  azimuth 
circle  observations.  Corrections  because  of  the  variability 


i 


of  the  measure  of  time  were  not  applied  in  any  case.  As  in 
the  case  of  Mars,  the  &S  of  the  declinations  measured  on  the 
azimuth  circle  are  different  from  the  values  published  in  the 
catalog,  because  the  latter  are  already  corrected  due  to  the 
systems  corrections  which  are  only  derived  here. 


A division  of  the  material  by  the  various  values  of  X.  -*c  gj  was 
not  made  in  the  case  of  Jupiter  and  Saturn  because  of  the  in- 
significance of  the  influence  of  the  corrections  of  the  earth 
orbit.  All  observations  which  were  made  during  an  opposition 
were  summarized  in  a normal  locus.  Tables  3 and  U present  the 
normal  loci  thus  obtained. 

For  the  determination  of  the  unknowns,  the  approximation  method 
described  for  the  Mars  observations  was  used,  consisting  of  an 
iteration  of  teh  equations  (36)  and  (37).  In  the  case  of 
Jupiter,  finally,  the  right  ascension  of  the  year  19^8  had  to 
be  omitted  because  it  could  not  be  reconciled  with  the  other 
observations  in  any  manner.  The  specific  observations  were 
not  made  with  the  meridian  circle,  but  with  a small  passage  in- 
strument, where  large  errors  are  possible  in  individual  cases, 
tg  = 19*H.O  and  10  years  as  the  time  unit  was  selected  in  the 
terms  with  a^.  Both  in  the  case  of  Jupiter  and  Saturn,  the 
coefficient  b^  resulted  only  with  considerable  uncertainty. 

The  final  numerical  results  were 

Et  a,  al  a,  6,  bt  6,  Du  Dc 

Jupiter  —2:02  +o"73  +o’.'55  — o'.'53  — <>''29  +o"2i  — 0:33  +o.'o5  — o"66 
Saturn  — 1.78  +2.37  +0.13  +0.29  ^+3.14  +0.03  — 0.25  +0.02  — 1.17 

The  evaluation  of  these  numbers  will  be  carried  out  in  the  next 
section  in  connection  with  a discussion  of  the  results  of  all 

planets . 


Table  3.  Normal  Loci  for  Jupiter 
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Table  U.  Normal  Loci  for  Saturn 


Rektaszensioneir  01 


da 


Deklinationen  •Jy' 


1.  Right  ascensions 

2.  Declinations 

3.  Time 

U.  Location 


Because  of  the  short  orbital  arc,  a determination  of  the  ele- 
ments for  Uranus  and  Neptune  was  impossible.  Only  the  decli- 
nations were  illustrated  through  a power  series  in  accordance 
with  time 

ja  = a + j(<_,o)  + e for  Munich  observations  ) (39) 

J<5  = a +-  b (t  — 1„)  4-  c (< — <o)s  + ^ for  Canberra  observations  ) 

Here  again,  Tq  = 19^1.0  and  10  years  was  selected  as  the  time 
unit.  Following  some  experiments  with  Neptune,  the  square  co- 
efficient c proved  to  be  indeterminable  and  was  then  set  equal 
to  zero.  The  term  d,  which  occurs  in  the  observations  in  Can- 
berra, is  equal  to  the  difference  D^  - Dc  of  the  two  equator 
corrections.  In  the  same  manner  as  for  Jupiter  and  Saturn,  all 
observations  made  during  an  opposition  were  summarized  to  a 


normal  locus; 


~2U 

these  normal  loci 
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are  shown  in  Table  5. 
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Key : 

1.  Time 

2.  Location 


These  values  were  subjected  to  an  adjustment  in  accordance  with 
formulas  (39):  the  following  results  were  obtained 

Uranus  NYptun 

a 1 '.'2.S  i i} ”42  -r  i Ji’ri:0  !7 

b — 0.27  i 2.43  i 0.25 

r — 0.33  ± 1.45  0.00  (anci-noimnin)  (assumed) 

d ± 0.75  - 0.41  ~ 0.32. 

The  values  d provide  additional  information  about  equator  cor- 
rections, though  very  unreliable,  and  were  used  in  the  catalog 
of  the  azimuth  circle  observations  (6)  in  this  manner. 

8.  DISCUSSION  OF  THE  RESULTS  OF  ALL  PLANETS  AND  COMPARISON 
WITH  THE  RESULTS  OF  THE  COMPUTER 

The  corrections  of  all  orbit  elements  and  the  system  improve- 
ments D^ , D^  and  E must  be  determined  from  the  values  of  the 
unknowns,  which  were  obtained  for  the  individual  planets.  The 
results  found  by  Petri  in  the  subsequent  work  concerning  the 
observations  of  the  small  planets  Ceres,  Pallas,  Juno,  and  Vesta 
can  also  be  used  for  this  calculation. 


The  following  values  have  resulted  for  the  equator  corrections: 


Planet 

* Munchen  (DM) 

Canberra  (Dc) 

^Munich 

Mars 

+o'-02 

+o"3i 

Jupiter 

+ 0.05 

-—0.66 

Saturn 

-r  0.02 

—1.17 

Ceres 

—>>.32 

—0.91  (1954) 
—0  37  (1055) 

**  undet 

Pallas 

unbest. 

O.JO 

Juno 

-fo.09 

4-0.73 

Vesta 

— 0.21 

— 051 

The  amounts  for  the  equator  correction  for  Munich  are  small 
and  scatter  little;  those  for  Canberra  are  large  and  scatter 
considerably.  It  is  shown  that  D^.  is  largely  a function  of 
tie  zenith  distance;  for  this  reason  two  separate  Dc  values 
were  derived  for  Ceres  , which  was  observed  in  Canberra  as  the 
only  planet  in  two  oppositions  with  significantly  different 
declination.  In  the  catalog  of  the  azimuth  circle  observation 
(6),  a correction  was  derived,  a.o.  from  the  results  of  the 
planet  observations  because  of  the  higher  terms  of  the  law  of 
curvature . 

The  determination  of  E was  faced  with  the  difficulty  that  the 
coefficient  b^  could  not  be  determined  in  most  cases  and  was 
only  unreliable  in  any  case.  For  this  reason,  the  second  equa 
tion  (31*)  was  utilized.  For  the  four  unknowns  A J,  AK , Ai  and 
E,  the  equations 

Pi  AJ  + AK  + £ de  = a, , 

qidj — p\dK — tj  de  = 6,, 

?i  dj  — pi  dK  — r,  de  + E tg  e = b,  — E0  tg  e 

are  thus  available,  which  must  be  established  for  each  group, 
into  which  the  material  of  the  observations  was  subdivided  c or 
responding  to  the  value  of  ■*--*  If  the  value  of  b2  is  un- 

certain or  not  sufficiently  reliable,  the  second  equation  is 
omitted  for  the  particular  group. 

In  the  case  of  all  planets,  the  resolution  in  accordance  with 
all  four  unknowns  was  so  significantly  uncertain  that  the  in- 
direct approach  was  used.  Under  the  first  hypothesis  £e.=  E = 
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temporary  corrections  A J and  AK  were  determined  for  each  planet 
and  the  remainders  were  again  solved  for  the  unknowns  and  E. 
Because  of  too  great  an  uncertainty  of  the  coefficients,  this 
calculation  was  not  performed  for  Juno.  For  Juprt-er  and  Saturn, 
E was  calculated  from  the  formula 

E tgf  - A,  — /-„  4:  f — I't 

which  is  evident  in  accordance  with  (Uo).  Values  of  Af  cannot 
result  in  the  case  of  Jupiter  and  Saturn  because  the  corrections 
of  the  earth  orbit  were  here  neglected  from  the  beginning.  The 
following  results  were  obtained: 

Table  6.  Results  for  4r  and  S 

I'lanct:  Jupitur  S.iturn  Mars  i rr«s  j'nll.js  \r^id  Mittrl  \Jean 

« » i - * * 

.if  — -f 0.03  —0.19  ! ! o-m  40.17  40 11*0.13 

£ | 4 o’tiS  J1T13  — 0.01  —0.12  I 40.4 1 40.73  4o.34io.il 

Because  of  the  lesser  accuracy  of  the  observations,  Jupiter 
and  Saturn  were  assigned  half  the  weighting  in  the  determina- 
tion of  the  mean  values  last  shown  in  Table  6.  The  two  mean 
errors  are  undoubtedly  too  optimistic  because,  in  the  method 
of  approximation  which  is  used,  the  values  of  the  unknowns  be- 
come quite  accurate,  but  the  mean  errors  become  too  snail. 

With  the  values  of  and  E,  which  were  obtained,  equations 
(Uo)  were  again  solved  for  the  two  unknowns  AJ  and  A K for  each 
planet. 

The  indirect  method  also  had  to  be  applied  for  the  determina- 
tion of  the  other  orbit  elements  in  accordance  with  formulas 
(32)  by  first  deriving  preliminary  values  of  the  planet  ele- 
ments for  each  planet,  while  ignoring  the  correction  quantities 
of  the  earth  orbit  and  then  subsequently  determining  the  earth 
orbit  quantities  from  the  remainders.  The  values 

Ah'  = -+-o"o7  i o'.'oS  , 

Ah'  --  +0.01  ± 0.05  , 

resulted,  of  which  the  mean  errors  simulate  a degree  of  accu- 
racy whi c h is  too  optimistic  for  reasons  which  have  already 
been  discussed. 
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In  the  determination  of  AL',  the  prerequisite  of  the  approma- 
tion.  method  is  not  fulfilled,  according  to  which  the  unknown 
AL',-  which  is  neglected  in  the  first  approximation,  must  he 
small  with  respect  to  the  others.  For  this  reson,  the  first 
two  equations  (32)  were  set  up  for  each  planet  and  Alq  and 
were  eliminated  from  the  normal  equations;  the  following  four 
final  equations  resulted  for  AL' 


Mars 

+0.17  AL' 

= ~r°"33 . 

Ceres 

+°P3 

= -f  O.IO  , 

Pallas 

+0.01 

= +0.07 . 

Vesta 

+0.05 

= +0.03 . 

Through  the  addition  of  these  four  equations,  the  result  was 
A L'  = +2.0U  ± 1"68 . The  corrections  of  the  earth  orbit,  which 
were  found,  were  inserted  in  the  original  equations  and,  from 
them,  in  a second  approximation,  the  final  corrections  of  the 
orbit  elements  of  the  planets  were  determined. 

The  definite  results  are  cited  in  the  subsequent  work  of  Petri 
in  Table  10,  which  also  gives  the  mean  errors.  The  same  table 
provides  the  results  of  a fitting  of  the  entire  material  of  the 
observations  by  all  unknowns  in  one  fell  swoop,  obtained  with 
the  PERM  computer  of  the  University  Munich;  these  results  were 
made  available  to  us.  In  general,  the  comparison  shows  good 
agreement.  Only  in  three  cases  (ALq  for  Pallas,  A u-  and  Aq;  for 
Saturn),  the  deviation  is  greater  than  the  mean  error. 

However,  the  mean  errors  determined  by  the  computer  are  so 
great  that  the  values  of  the  unknowns  cannot  be  guaranteed  in 
almost  all  cases.  It  could  be  expected  from  the  beginning  that 
material  consisting  of  somewhat  more  than  500  individual  obser- 
vations cannot  do  justice  to  the  complete  determination  of  all 
unknowns.  The  large  mean  errors  are  significantly  a consequence 
of  this  excess  demand  on  the  material  of  the  observations. 


Yet,  the  calculation  carried  out  with  the  PERM  computer  has 
the  considerable  value  of  having  proved  the  principal  solubi- 
lity. of  as  extensive  a task  as  the  simultaneous  determination 
of  50  unknowns.  In  contrast,  the  approximation  method  describe 
in  this  paper  proves  the  possibility  that,  even  without  the 
use  of  a computer,  which  is  often  not  available  or  can  be  used 
only  with  the  use  of  relatively  significant  resources,  the 
problem  can  be  solved  manually,  when  the  requirements  for 
accuracy  are  moderate. 

The  uncertainty,  which  is  inherent  in  the  calculation,  in  the 
case  of  the  equinox  correction  E,  becomes  particularly  clear. 
Because,  on  the  one  hand,  only  the  quantity  E tanfc-can  be 
determined  with  complete  accuracy  and,  on  the  other  hand,  in 
the  material  of  the  observations,  the  right  ascensions  are 
small  in  number  and  have  an  unfortunate  distribution,  the  sig- 
nificant uncertainty  of  E could  be  expected.  The  discrepancy 
between  the  formal  value  +1"30  - 1"36,  found  with  the  PERM  and 
the  value  +0"3i*  - 0"21  is  so  great  that  neither  one  of  them 
can  lay  claim  to  objective  accuracy.  nevertheless,  the  author 
believes  that  the  second  value  E = +0"3^  is  reliable;  it  has 
good  agreement  with  other  determinations  of  E and  it  can  claim 
the  fact  of  experience  that  the  method  of  approximation  which 
is  used  mostly  results  in  reasonably  correct  values  for  the 
unknowns,  but  has  mean  errors  which  are  too  small. 

The  fact  that  the  equator  corrections  and  Dc  which  were 
found  have  resulted  in  a good  agreement  of  the  declination 
systems  of  the  Munich  azimuth  circle  for  the  observations  in 
Munich  and  Canberra  is  a very  powerful  argument  for  the  usabi- 
lity of  the  described  approximation  method.  Thus,  it  is  shown 
that  the  system  (31)  is  suitable  for  a simple  evaluation  of 
meridian  observations  of  the  outer  planets,  which  is  sufficient 
ly  accurate  with  a limited  amount  of  material. 


! 


-29- 

BIBLIOGRAPHY 

1.  Cleraence,  G.  M.:  The  Value  of  Minor  Planets  in  Meridian 
Astronomy,  ASTRON.  JOURNAL,  51* : 10(  19^8  ) . 

2.  Annals  of  the  Cape  Observatory,  Vol . 8,  Parts  k (1915) 
and  5 (1921). 

3.  Numerow,  B.:  On  Providing  a Fundamental  System  of  Weaker 
Stars,  ASTRON.  NACHR. , 260:305(1936). 

4.  Bauschinger,  J.:  The  Orbit  Determination  of  Celestial 
Bodies,  2nd  edition,  p U50(1928). 

5.  Newcomb,  S.:  Formulae  and  Tables  for  Expressing  Corrections 
to  the  Geocentril  Place  of  a Planet  In  Terms  of  Symbolic 
Corrections  to  the  Elements  of  the  Orbits  of  the  Earth  and 
Planet,  ASTRON  PAPERS,  Vol.  II,  Part  I (1891). 

6.  Schmeidler,  F.:  Measurements  of  Fundamental  Declinations 
in  Both  Hemispheres,  VEROEFF.  STERNW.  MUENCHEN,  Vol.  k. 

No.  22  (1957). 

7.  Labitzke,  P. : Observations  of  Planets  on  the  six-inch 
Meridian  Circle  of  the  Munich  Observatory,  ASTRON.  NACHR.. 
275  : 133  ( 19**7  ) - Published  in  STERNW.  MUENCHEN,  Vol.  3,  Ho  . 5 . 


8.  Labitzke,  P. : Right  Ascension  Observations  of  Planets  on 
the  Broken  Transit,  ASTRON.  NACHR.,  280:259(1952)  - 
Published  STERNW.  MUENCHEN  Vol.  U,  Ho.  9. 


